PERIODIC SOLUTIONS OF PERIODICALLY PERTURBED PLANAR 
AUTONOMOUS SYSTEMS: A TOPOLOGICAL APPROACH t 



Mikhail Kamenskii 1 , Oleg Makarenkov 1 and Paolo Nistri 2 

department of Mathematics, 
■ Voronezh State University, Voronezh, Russia. 

, E-mail: Mikhail@kam.vsu.ru, E-mail: omakarenkov@kma.vsu.ru 

2 Dipartimento di Ingegneria dell'Informazione, 
Universita di Siena, 53100 Siena, Italy. 
jy^ • E-mail: pnistri@dii.unisi.it 

OO 

Abstract 

' Aim of this paper is to investigate the existence of periodic solutions of a nonlinear planar autonomous system having a limit cycle 

xq of least period To > when it is perturbed by a small parameter, Tj — periodic, perturbation. In the case when Tq/Ti is a 
^ (— I rational number l/k, with I, k prime numbers, we provide conditions to guarantee, for the parameter perturbation e > sufficiently 

small, the existence of HTq— periodic solutions x e of the perturbed system which converge to the trajectory xq of the limit cycle 
as e — ► 0. Moreover, we state conditions under which T = klTo is the least period of the periodic solutions x £ . We also suggest a 
simple criterion which ensures that these conditions arc verified. Finally, in the case when To /Ti is an irrational number we show 
the nonexistence, whenever T > and e > 0, of T— periodic solutions x e of the perturbed system converging to xq. The employed 
i methods are based on the topological degree theory. 
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1 Introduction 

The paper is devoted to the study of the existence and of the behaviour as e —* of periodic solutions 
of a perturbed autonomous differential system in R 2 of the form 

x = ip(x) +e<j>{t,x), (1) 

where tp : R 2 — > R 2 is a twice continuously differentiable function, <fi : R x R 2 — > R 2 is continuous and 
Ti-periodic with respect to the first variable and e > is a small parameter. 
To be specific, we assume that at e = the autonomous system 

x = ^(x), (2) 

has a limit cycle xq = xo(t), t > 0, of least period Tq satisfying the following condition 
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methods", by GNAMPA-CNR, by RFBR grants 02-01-00189, 05-01-00100, by U.S.CRDF - RF Ministry of Education 
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(Aq)- the linear system 



y = ^(x (t))y (3) 
does not have 2To-periodic solutions linearly independent with xo(t). 

Furthermore, we assume 

T I 

(Ai)- there exist I,fc£N such that — = — , with I and k prime numbers. 
This paper addresses the following problems: 

1. To provide conditions on the function (p which guarantee the existence of > such that system 
(1) has a T '-periodic solution x £ , with T = T(Tq,Ti) and e G (0, £*) satisfying the property 

x e (t) — ► xq as e — > 0, whenever t E [0, T], 

where xq = {x G R 2 : x = a;o(i), i G [0, To]} is i/te trajectory of the limit cycle of (2). 

2. To find an explicit estimation of > 0. 

3. To investigate the existence of periodic solutions for system (1) in the case when condition (A±) is 

iji 

not verified, i.e. when — is an irrational number. 

J-i 

Several contributions have been made toward solving the problem of the existence and the nature 
of periodic solutions of a periodically perturbed autonomous system which are close to the periodic 
solution of the autonomous system. We refer to the papers [20], [21], [19] and to the memoir [22] for 
first and second order perturbed autonomous systems respectively. The employed methods are based 
on perturbation and implicit function techniques which require more regularity on the function 4> of 
that required in this paper. In [20], under assumptions similar to (Aq) and (Ai) the author provides 
sufficient conditions for the existence and stability of periodic solutions of a perturbed autonomous 
system in R n . We would like to point out that in this paper we are interested in periodic solutions 
of general form, that is without any condition on their rate of convergence to xq with respect to e. 
In such a generality the uniqueness of the periodic solutions is not guaranteed as in [20], where a 
special class of periodic solutions, converging to xq at the same rate of e, is considered. Similar results 
for the existence of D— periodic functions, namely for functions whose derivative is periodic, can be 
found in [9] and [10]. Furthermore, in [3], [4] and [5] analogous existence and stability results are 
established under different assumptions. Specifically, in [3], under the assumption that 1 is the only 
characteristic multiplier of the linearized system, which is a q-th root of unity, for q G N, it is proved 
that for any p G N there exists a parametrized family of periods and corresponding periodic solutions 
of the perturbed autonomous system such that at e = the period is given by qT^/p. In [4] and [5] 
the behaviour, with respect to the period of the disturbance, of the periodic surface described by a 
functional equation and associated to the perturbed autonomous system is studied. Finally, in [11] 
and [12] existence and multiplicity results of periodic solutions of a periodically perturbed autonomous 
second order equation, defined on a manifold, were proved by means of topological methods. An 
example of periodically perturbed autonomous systems arising in Mathematical Biology is given in 
[13]- 

An other approach to the study of the existence of periodic solutions of perturbed autonomous systems 
based on topological degree and index theory is presented in [8] and [1] respectively. In these papers 
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the perturbation is not necessarily small, in the sense that it does not contain small perturbation 
parameters, the period T is fixed and the autonomous system may have periodic orbits of period less 
than T. 

In this paper to solve the proposed problems we use a different approach which was previously intro- 
duced and successfully employed by the authors in [14], [15] and [16] to treat the existence problem 
of periodic solutions of perturbed nonautonomous systems of prescribed period. Specifically, in [14] 
we consider the nonautonomous system of differential equations described by 

x = ip(t, x) + e<f>(t, x) , (NA) 

where 4>, tp '■ R x R n — > R n are continuously differentiable, T-periodic with respect to time t, functions 
and e is a small positive parameter. 

The proposed approach is based on the linearized system associated to (NA) 

y = ^(t, n(t, o, 0)y + <P(t, n(t, o, 0), (lna) 

where £ G R ra and J7(-,io,£) denotes the solution of (NA) at e = satisfying x(to) = £. Specifically, 
consider the change of variable 

Z (t) = n(o,t,x(t)), (cv) 

and the solution rj(-, s, £) of (LNA) such that y(s) = 0. If there exists a bounded open set U C R n such 
that n(T, 0, £) = £ for any £ G dU, and r){T, s, f ) - r?(0, s, £) / 0, for any s G [0, T], and any £ G dU. 
Then (NA) has a T-periodic solution for e > sufficiently small provided that T(rj(T,0,-),U) / 0. 
Here r(i ? , U) denotes the rotation number of a continuous map F : U — > R 2 . Observe that T(F, C7) 
coincides with deg(F, f/, 0), the topological degree of F at relative to U. In what follows, we will 
omit the point in the notation for the topological degree. 

The advantage of the proposed approach as compared with the classical averaging method, which is 
one of the most useful tools for treating the existence problem of periodic solutions of nonautonomous 
periodic systems, mainly consists in the fact that in order to use this second method for establishing 
the existence of periodic solutions in perturbed systems of the form (NA) one must assume that the 
change of variable (CV) is T-periodic with respect to t for every T-periodic function x such that 
ri(0, t, x(t)) G U, for any t G [0, T], instead that only on the boundary of the bounded open set U. 
The same assumption is necessary in vibrational control problems, [2] and [23], to reduce the considered 
system to the standard form for applying the averaging method. For an extensive list of references on 
this topic see [7]. 

Our approach has been also employed in [15] to prove the existence of periodic solutions for a class of 
first order singularly perturbed differential systems. 

Furthermore, in [16] we have considered two nonlinear small periodic perturbations of ip(t, x) with 
multiplicative different powers of e > and we have proved the existence of T— periodic solutions of 
the resulting system. We have also showed the presence of the so-called frequency pulling phenomenon 
in the case of a special class of planar autonomous systems when they are perturbed by periodic terms 
with period close to that of the periodic solution of the autonomous system, whose existence is assumed. 
The aim of this paper is at solving problems 1-3 by using the approach outlined above, extending in 
this way its application to the study of the existence of periodic solutions and of their behaviour for 
general, periodically perturbed, planar autonomous systems having a limit cycle. In particular, we 
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are interested in investigating the relationship between the period of the limit cycle and that of the 
nonautonomous perturbations in order to have periodic solutions of (2) with the period expressed in 
terms of the previous ones. 

To the best knowledge of the authors the present paper treats for the first time the problem of 
evaluating > 0, i.e. problem 2. stated before, for the application of the averaging principle via 
topological methods. 

The paper is organized as follows. In section 2 we show our main result: theorem 1 which provides 
sufficient conditions to ensure, for e > sufficiently small, the existence of MTq— periodic solutions x £ 
to (2) converging to xq as e — > 0. In theorem 2 we show, under a mild extra assumption, that klT is 
the least period. 

In section 3 we present a simple method to ensure that the conditions of theorem 1 are satisfied. 
Finally, in Section 4 we establish a nonexistence result: theorem 4, in the case when Tq/T\ is irrational. 
Precisely, we show that, for any T > and e > 0, there is not a T— periodic solution x £ of (2) such 
that x £ — > xo as e — ► 0. An example illustrating the existence result is also provided. 

2 The main result 

Throught this section we assume that T = HTq and we will denote by FL the derivative of the 
function F with respect to the i-th argument. Let x{t) = Q(t, 0,£) be the solution of system (2) 
satisfying x(0) = £. Consider the following auxiliary system of linear ordinary differential equations 

y = ^'(n(t,o,0)y + (t>(tMt,o,0)- (4) 

and let y(t) = rj(t, s, £) the solution of (4) satisfying y(s) = 0. Observe that if ip = then 

T 

V (T,s,0-v(0,s,0= J Hr,0dr. 

o 

Therefore, in this case, the function — (rj(T,s,£) — rj(0, i s the average on the interval [0, T] of 
the function <p with respect to the first variable. This function is the basis of the classical averaging 
method, one of the most relevant tools to investigate the existence of periodic solution of (NA) when 
t/j = (see [6]). 

We can prove the following preliminary result. 
Lemma 1 

t 

s 

where 

^(t,o = %)(o,tMt,o,omMt,o,0)- (5) 

Proof: Observe that the matrix tt'^(t, 0,£) is the fundamental matrix, satisfying J7'^(0, 0,£) = /, 
for the linear system 

y = VW,0,0)y. 
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Furthermore, (il'^{t, 0,£)) = ^( 3 )(0, t, Q(t, 0, in fact, by deriving with respect to £ the identity 

n(0, t, f2(i, 0, f )) = f , whenever £ G R 2 , 

we obtain 

n[ 3) (0,t,n(t,0,0)n[ 3) (t,0,0 = I, whenever £e R 2 . (6) 
Therefore, by the variation of constants formula for the nonhomogeneous system (4) we obtain 

t 

s 

t 

s 

□ 

Since the trajectory xq is a Jordan curve, its interior U is a bounded, simply connected, open set of 
R 2 . In order to prove, by the proposed approach, the existence of periodic solutions to (1) converging 
to x as e — > 0, we introduce a family of open sets W 1 (U) as follows 



W,(U) 



U\(dU+\i\B) if 7 <0, 
U U (dU + |7|B) if 7 >0, 



where B C R 2 is the open unit ball, thus we have 

B 7 (dU) := ((W 7 (U) U U)\W 7 (U)) (J ((W 7 (E0 U f/)\[7) -► x as 7 0, (7) 
where <9?7 = xo- 

Definition 1 Define the following positive constants 

M 1 = max ||$(i,£)||, 

te[o,T], £eB 7 (dU) 

M!y = max ||$(2)(U)I|, 

te[o,T], feB 7 (at/) 

L; = max (T, 0,0||, 

ZeB^(dU) 

L>> = max_||^ 3)(3) (T,0,e)||, 

£€B~/(dU) 

Moreover, let 70 > suc/t i/mi system (2) does not have To-periodic solutions with initial condition 
belonging to the open set B_ l0 (dU) [j B^ Q (dU). 
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Observe that condition (Ao) guarantees the existence of the constant 70 in the previous definition. 
We can now formulate the main result of the paper. 

Theorem 1 Assume (Ao), (A\) and 

(A 2 )- f](T, s, f ) - f](0, s, f ) / 0, for any s G [0, T], and any £ G dU, 

(A 3 )- deg(^(0,T,-),l7)^l. 
Then for every —70 < 7 < 70, 7^0 and for 

\ K o } 

< e < min < t r-, — — — — > := e 7 (8) 

\ T2M 7 (M> + x 2.\/- /.» + .!/.'/.'_ ) ™ 7 (1 j 
system (1) has a T -periodic solution x £ belonging to the set 

[x e R 2 : x = fi(t,0,0, * G [0,21, C G S t(^)} • 
To prove this theorem we need the following lemma. 

Lemma 2 Assume (Ao). Then there exists 71 > such that for —71 < 7 < 71 and 7 / we have 

deg(I-n(T,0,-),W 7 (U)) = 1. 

Proof: Condition (Ao) ensures that the characteristic multiplier of Q'^(T, 0, x(0)) different from +1 
is not equal to —1. Therefore the limit cycle xq is either asymptotically stable or unstable, thus in a 
sufficiently small neighborhood of xo the map £ — > fJ(T, 0, £) does not have fixed points different from 
those belonging to xo- in particular, 

^n(T,0,0 for any £ G cW 7 (£Z) (9) 

and I7I / sufficiently small. Hence by Corollary 2 of [8] we have that 

deg(J - n(T, 0, •), W 7 (U)) = deg(V, W 7 (U)), 

On the other hand from (9) we have that 

deg(^, W 7 (C/)) = deg(V>,£/) 

and since the vector field ip is tangent to at any point, we get (see, for example, Theorem 2.3 of 
[18]) 

degty,U) = 1. 

□ 

The following result is a direct consequence of Lemma 2. 
Corollary 1 Under condition (Ao) we have that 

deg(7 - fi(T, 0, •), W 7 (17)) = 1, for every - 7o < 7 < 7o, 7^0. (10) 
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Proof: First we prove (10) for < 7 < 70. By lemma 2 there exists < 7* < 70 such that 
deg(I — £l(T, 0, •), W 7 „ (U)) = 1. From definition 1 it follows that the constant 70 has the property that 
the vector field I — fi(T, 0, •) is not degenerated on the boundary of the set W 7 (C/) for any < 7 < 70. 
Therefore, for any < 7 < 70, we have 

deg(J - n(T, 0, •), W 7 (U)) = deg(J - Q(T, 0, •), W 7 , (U)) = 1. 

The same arguments apply for —70 < 7 < 0. 

□ 

We are now in the position to prove theorem 1. 

Proof: Denote by C([0,T],R 2 ) the Banach space of all the continuous functions defined on the 
interval [0,T] with values in R 2 , equipped with the sup-norm. Consider in (1) the change of variable 

x(t) = n(t,o,z(t)). (11) 

For every z G C([0, T],R 2 ), (11) defines uniquely x G C([0, T],R 2 ) with inverse given by 

z(t) = n(0,t,x{t)), t G [0,T]. (12) 

Therefore, the function x is the solution of the system (1) if and only if the function z defined by (12) 
satisfies the differential equation 

n' (1) (t, 0, z(t)) + n{ 3) (t, 0, z(t)) z{t) = £ <f>(t, n(t, o, z{t))) + o, z(t))). (13) 

By the definition of 0, z(t)) we have 

f2' (1) (*,0,z(t)) =^(*,0,z(t))). (14) 

Moreover, by using (6) and (14) we can rewrite system (13) in the following form 

z(t)=e^(t,z(t)), (15) 

where $ is defined by (5). Observe that (11) and (12) define a homeomorphism between the solutions 
of systems (1) and (15). 

Consider an arbitrary T-periodic solution x of system (1), we have 

z(0) = O(0, 0, x(0)) = s(0) = x{T) = Q(T, 0, z{T)). 

Therefore the problem of the existence of T-periodic solutions to system (1) is equivalent to the 
problem of the existence of zeros of the compact vector field G £ : C([0, T], R 2 ) — ► C([0, T], R 2 ) defined 
by 

t 

G £ (z)(t) = z(t) - n(T, 0, z(T)) - e J $(r, z(r))dT, t G [0, T]. 

o 

Define the set 

Z = {z G C([0,T],R 2 ) : z(t) G B^dU) for any t G [0,T]} . (16) 
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Consider the auxiliary compact vector field 

G M = I-A £ : C([0,T],R 2 ) - C([0,T],R 2 ), 

where 

T 

A £ (z)(t) = fi(T, 0, z(T)) - e J $(r, z(r))dr, t G [0, T], 

o 

Let us show that for any e > satisfying (8) the compact vector fields G £ and G\ )£ are homotopic on the 
boundary of the set Z. For this, define the following homotopy F £ : [0, 1] x C([0, T], R 2 ) — ► C([0, T], R 2 ) 
joining the vector fields G £ and G\ j£ : 

a{X,t) 

F £ (X,z)(t)=z(t)-n(TAz(T))-e J $(r, z{r))dT, t £ [0, T], 

o 

where a(A, t) = At + (1 — A)T. Let us show that for any e > satisfying (8) the homotopy F £ does not 
vanish on the boundary of the set Z. Assume the contrary, thus for some e > satisfying (8) there 
exist z £ G dZ and A e 6 [0, 1] such that 



a(A e ,t) 



; {t) = n{T,0,z £ (T))+e J Q(T,z e (T))dT, t e [0,T]. 



(17) 



From the fact that z £ G it follows the existence of t e G [0, T] such that .z e (£ £ ) G d{B 1 {dU)) . By 
the definition of the set B^{dU) either 



or 



Z £ {te) g dw 7 (u), 

Ze(t £ ) G 8U. 



(18) 

(19) 



By using (8) we have the following estimate 



\z £ (t £ )-n(T,0,z £ (t £ ))\\ = 

\z e {te) - Sl(T, 0, z £ {T)) + n(T, 0, z £ {T)) - (l(T, 0, * e (t e ))|| 

a(\ e ,t e ) 

J $(r, z £ (r))d T + Q(T, 0, z £ (T)) - fi(T, 0, « e (t e )) 
o 

< eTM 7 + eL' 1 TM 1 < K 7 , 



< 



(20) 



which contradicts the definition of the constant i^T 7 in the case when (18) holds true. Thus (18) cannot 
occur. Assume now (19), by lemma 1 and the fact that, in this case, z £ (t £ ) = Q(T,0,z E (t £ )) we have 

V (T,a(\ £ ,t £ ),z £ (t £ ))- V (0,a(X £ ,t £ ),z £ (t £ )) = 



(n{ 3) (T,0,Z £ (t £ ))-l) f <S>(T,Z £ (t £ ))dT+ U(T,Z £ (t £ ))dT 



a(\ £ ,t e ) 
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z £ (T) - z e {te) Sl(T, 0, z £ {T)) - n(T, 0, z £ (t e )) 



1 

-J §(T,z e (r))dT- 



1 1 

(l-n{ 3) (T,0,Z £ (t £ ))) J <S>(T,Z £ (t £ ))dT + J <S>(T,Z £ (t £ ))d7 

[ ' £i )] e \ [n{ 3) (T,0,^eKz £ (T)-z £ (t £ ))] 2 ) 



Ot(\ £ ,t € ) 



-(l-n{ 3) (T,0,Z £ (t £ ))) J <S>(T,Z £ (t £ ))dT + J<Z>(T,Z £ (t £ ))dT- 



a(\ £ ,t e 
a(\e,t e ) 





a(X e ,te) 



J <Z>(T,Z £ (T))dT = J <Z>(T,Z £ {t £ ))dT - J <Z>(T,Z £ (T))d7 



' [jy (3) (T,(U l£ ) J *(r,z e (r))dr]i X 

a(\ e ,t s ) 

[jy (3) (T,0,k £ ) J 0>(r,z e (r))dT] 2 

a(A £ ,t £ ) 



+ 



+ 



1 [n[ 3) (T,0,z £ (t £ )) J $(r,z e (r))dr]i N 

a(Ae,t e ) 



[n' (3) (r,0,z e (t e )) / $(r,z £ (r))dr] 2 



a(A £ ,t s ) 
T 



/ 



-O' (3) (r,0,^(t £ )) | $(T,Z £ (T))dT+ n{ 3) (T,0,Z £ {t £ )) J $(T,Z £ (t £ ))dT, 



(21) 



a(A e ,t £ 



«(A £ ,t £ ) 



where ^ e G [z £ (t £ ), z £ (T)\i and i = 1, 2, denotes the i-ifo component of the vector v. By using the 
following estimates 



\z e {t e ) - < \\z £ {t £ ) - z £ {T)\\ < 



$(r,z e (r))dr 



a(\ e ,te) 



< eTM~ and 



max ||$(t,z e (t e )) - $(t,z e (t))|| < M^max \\z £ (t £ ) - z £ (t)\\ < eTM y M^ 

together with (8) we obtain from (21) 

\\ V (T, a{X £ ,t £ ), z £ {t £ )) - 7/(0, a{X £ ,t £ ),z £ {t £ ))\\ < 
< eT 2 M^M^ + eV2 T 2 (Mj) 2 L" + eT 2 M 7 M 7 L 7 = 

= eT 2 M 7 (M 7 + v / 2M 7 L^ + M 7 L 7 ) < K , (22) 

which contradicts the definition of the constant Kq. Hence both (18) and (19) cannot occur and so for 
e G (0, e 7 ) the homotopy F £ does not vanish on the boundary of the set Z and so the compact vector 
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fields G £ and G\ )£ are nomotopic on the boundary of the set Z. In particular, 

G h£ (z) + 0, for any z G 8 (zf| C const ([0, T], R 2 )) C 8Z, and any e G (0,e 7 ), (23) 

where C coris t([0, T], R 2 ) denotes the subspace of the space C([0,T],R 2 ) consisting of all the constant 
functions defined on the interval [0, T] with values in R 2 and d [Z f] C conS (([0, T], R 2 )) is the relative 
boundary of the set Zf] C cons t([0, T], R 2 ) with respect to the subspace C con <rf([0, T], R 2 ). Therefore, 
by the reduction domain property of the topological degree, taking into account that A £ (dZ) C 
C const {[0,T],R 2 ), we obtain 

deg C ([o, n R2)(Gi ie ,Z) = deg Cconst([0iT]iR2) (g 1;£ , Zp|C const ([0,T],R 2 )) , (24) 

whenever e G (0,e 7 ). Observe, that z G Z |"| C cons i([0, T], R 2 ) is a solution of the equation G\ i£ z = if 
and only if £ = z is a solution of the equation Q £ £ = 0, where Q £ : R 2 — > R 2 is defined by 

T 

o 

Therefore, we have 



de gc const ([o,T],R 2 )(G'i, £ ,ZnC const ([0,r],R 2 )) = deg R2 (Q e ,.B 7 (0l7)) 
= sign( 7 ) (deg R2 (Q £ , W 7 (17)) - deg R2 (Q £ , U)) , for any e G (0,e 7 ) 
From (8), for any e G (0,e 7 ), we have 



(25) 



T 

<J>(r,£)dT 



o 



and so 

deg R2 (Q £ , W 7 (?7)) = deg R2 (Q , W 7 (f/)), whenever e G (0,e 7 ). 
By corollary 1 and the previous equality we obtain 

deg R2 (Q £ , Wj(U)) = 1, whenever e G (0, e 7 ). (26) 

Let us now calculate deg R2 (Q £ , U). For this, let <5i, £ : R 2 — > R 2 be defined by 

T 

QlA = ~ £ j *(r,S)dT, for any £ G Z f]C const ([0,T],R 2 ). 
o 

From (A2) we have that Q £ £ = Qi, £ £, for any £ G dU, since £ = fl(T, 0,£), and so for e G (0,e 7 ) we 
have 

deg R2 (Q £ ,£7) = deg R2 (Q li£ ,C/). (27) 
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Let us show that the compact vector fields Qi )£ and Qi,i are homotopic on the boundary of the set U 
for e G (0,e 7 ). For this, define the linear homotopy F^ e : [0, 1] x R 2 — > R 2 as follows 

T 

F M (A,0 = -(Ae+l-A) | $(r,Odr. 

o 

Arguing by contradiction, assume that for some AG [0,1], £ G <9?7 and e G (0,e 7 ) we have 

T 

(Ae + l-A) J $(r,OdT = 0, 
o 

and so 



T 

| $(r,|)dr = 0. 
o 

By lemma 1 this is a contradiction with condition (A2). Thus, again by lemma 1 



deg R2 (Q 1;£ ,C/) = deg R2 (Q 1;1 ,C/) =deg R2 (7?(0,r, ■),[/)■ (28) 
From (27) and (28) we obtain 

deg R2 (Q e , U) = deg R2 (r/(0,T, •), U), for any e G (0,e 7 ). (29) 

Finally, taking into account (24)-(26), (29) and condition (A 3 ) we get 

de gc([o,T],R2)(G e , Z) = deg C ([ 0)T])R2 )(Gi )£ ,Z) = 
= sign( 7 )(l-deg R2 (7 ? (0,T,.),£/))/0, 

for any e G (0, e 7 ). The solution property of the topological degree ends the proof. □ 

Remark 1. Observe that deg(ry(0, T, •), U) can be calculated by means of a well known formula (see 
e. g. [17], p. 6). In this case condition (As) takes the form 

n-l 

deg(rK0,T, ■),£/) = ^sign([r/(T,0,x (-))]i,^ + i,^)- 

•(signer, 0,x (# i+ i))] 2 - signfo(T,O,so(0i))] 2 ) / 1, 

where j = 1,2, denotes the j — th component of the vector field r] and 9i, i = 0, 1, n — 1, # n = #0j 
are the ordered roots of function [t/(0, T, xo(#))]i- 

Conditions for the existence of solutions to (1) in terms of the derivatives of the function [i](0, T, xo(-))]i 
at the points 9i, up to a suitable rotation, are given in [20]. 

We now prove the following result. 

Theorem 2 Assume that there exists ^ G xq such that T\ > is the least period of the function 
t — > </>(t,£). Assume condition (A\), if {x £ } ee ^ £ ^ are T— periodic solutions to (1) converging to xq as 

e — > 0, then MTq < T. 
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Proof: Assume the contrary, thus there exists T G (0, klT$) such that for every e > sufficiently 
small system (1) has a T-periodic solution x £ satisfying 

lim xJt) = xo(t + wo), 

for some wq G [0, To], thus T is a period of the function xo(t). On the other hand, Tq is the least period 
of Xo(t) and so we obtain T = uqTq for some no G N. Moreover, in (1) we have that the function 
t — > (p(t,x £ (t)) is T-periodic, since i(i) is T-periodic, and so for any i £ [0)^1 > we have 

<j>(t ,x £ {to)) = <t>(t + f , x £ (t + T)) = </>(t + T, x e (t )). 

Due to our assumption and the arbitrarity of to, by passing to the limit as e — > 0, we conclude that 
T = n\T\ for some ni G N. Therefore 

n T = mTi 

and so by {A\) there exists ri2 G N such that no = klri2, contradicting the fact that T < MTq. □ 



3 A method to verify assumptions (A2) and (A3) of theorem 1 

Although the topological degree (or rotation number) of the vector field in (A3) can be calculated by 
using the methods from [17], as recalled in remark 1, we provide here an alternative method to verify 
conditions (A2) and (A3) which turns out to be useful and of simple application. 
First, we introduce some notations. Recall that for the vector v G R 2 we denote by [v]i its i-th 
component, i = 1,2, for a G R 2 we put 

„-L ( -W2 ^\ „T _ ( i a h \ 



and for a, b G R 2 




Wi H 2 1 („h\-( ^1 ^! 

Mi M2 r {ab) ~{ [a} 2 [b} 2 



Moreover, y = y(t) will denote the solution of system (3) linearly independent with xo(t) and (s, 0) 
F(s,9) is the following function 

s 

F(s,6)= J (± (t) 2 /(r))- 1 ( /»(r-^xo(r))dT. 

s-T 



We can prove the following. 

Theorem 3 Let f : [0, To] —> R 2 be a function satisfying the following conditions 
(B 2 ) (F(s,e),f(9)) + 0, for any s G [0,T], 6 G [0,T ], 
(B 3 ) deg(iV, [0, T ]) + \, 

where N(0) = (y(0) T i (6 l ) ± ) f(0), G [0,T ]. Then assumptions (A 2 ) and (A3) of theorem 1 are 
satisfied. 
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To prove this theorem we need the following preliminary lemma. 
Lemma 3 We have that 

( V (T, s, x (9)) - 77(0, s, x (0)), N(0)) = (30) 

= *f (d(T,0)<f>(T-0,x (r)), < x o (0),N(0) > y(r) T + < y(9),N(9) > f (r) ± )dr, 
s-r+e v ' 



where 



Proof: Observe, that 

n' (3) (o, t, n(t, o, x (9))) = k(o, 6) (K(t, e))~\ 

where K(-, 9) is a matrix whose columns are linearly independent solutions of the differential system 

x = ip'(xo(t + 9))x. (31) 
Let us choose the matrix K(-,9) as follows 

K(t,9) = (x (t + 9) y{t + 9)). 

Thus 

n{ 3) (0,t,n(t,0,x (9))) = (x (9) y(9))(x (t + 9) y(t + 9))' 1 = 

= (x (9) y(9))( y{t ^ ± )d(r,9). 
\ x {t + 9) J 

Furthermore, we have 

(r)(T,s,x (9))-r)(0,s,x (9)),N(9)) = 

f2' (3) (0, r, fi(r, 0, x (9)))</>(t, x (t + 9))dr, N(9)\ = 

\s~T ' 



s-T 



J d(r.e) (x„(«) y(«)) y ^ T + + "s)-L ) *^^(r + D))dr.N(e^ 

<x o (0),iV| 
<y(9),N(6)> 



"s-T 
s 



= J Ur,9MT,x (T + 9)),(y(T + 9) T x (r + 9)^ I < )V 
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s 

= J (cI(t,0)<P(t,x o (t + 8)),< x o (6),N(0) >y ( T + 9) T +<y(9),N(8) > x q {t + 8) L ) dr 

s-T 

s+e 

= J (d(r, 0)<j,(r - 9, x (t)), < xo(6),N(e) > y(r) T + < y(9),N(9) > x (r) ± ) dr. 
s-T+e 

We can now prove Theorem 3. 

Proof: By lemma 3 and the fact that 

(x o (6),y(0) T ) = (y(6),i o (6) ± ), 

we have 

(r j (T,s,x (8))- V (0,s,x (8)),N(8)) = 
s+e 

J (d(r, O)0(r - 0, x (r)), [/(^ (i o (0), y(#) T ) y(r) T + 

+[/(0)] 2 (*„(*), xo W ± ) y(r) T + [/(e)}, (y(9), y(9) T ) ±0^ + 
+ [f(9)} 2 (y(9),x (9) ± )x (r) ± )dT = 



s+0 

J (d(r, O)0(r - 0, x (r)), (x (#), y(#) T ) ([/(l^ y(r) T + [f(9)] 2 ^(r)^) ) 
T+e 

s+e 

x (9) y(9) T ) J (d(r, 0)<f>(r - 9, x (r)), (y(r) T x (r) ± ) /(<?)) dr = 

s-T+6* 

i ((9) 2/(^) T ) / (d(r,0) ( ) #r - (9,x (r)),/((9))\ dr. 



□ 



Therefore, by the condition (B2) 

(r ? (T, S ,x o (0))-r / (O, S ,xo(0)),^)) / 0, for any s G [0, T], 8 € [0, T ], (32) 

that is the condition (A2) of Theorem 1 hold true. Condition (A3) of theorem 1 also follows from (32) 
by taking into account condition (B3) and lemma 1. □ 

T 

4 The case when — is irrational 

Ti 

In this section we assume the following condition 
(A[)- — is an irrational number. 
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We can prove the following result. 

Theorem 4 Assume that there exists £ G xq such that the function 4>(t,£) is not constant with respect 
to t. Assume (A[), then system (1) does not have T— periodic solutions x £ , e G (0,e), converging to 
xq as e — ► 0, whenever T > and £ > 0. 

Proof: Assume the contrary, thus there exists T > and £ > such that for every e 6 (0, £) system 
(1) has T-periodic solution x £ (t) satisfying 

limx e (t) = x (t + w ), 

for some wo G [0, To]. Hence T is a period for the function xo- On the other hand, since To is the 
least period of x$, we have that T = noTo for some no G N. Moreover in (1) we have that the function 
t — > (p(t,x £ (t)) is T-periodic, since x(£) is T-periodic, thus, for any to £ [0, T], we get 

^(io^eOo)) = 0(*o +pn T ,x £ (t +pn T )) = </>(i + T p ,x e (t )), whenever p G N, 
where T p = pnoT]( mod Ti). Condition (A^) implies that 

U^=[0,T 1 ] 

peN 

and so 

(j)(t,x £ (t )) = c £ , for any t £ [0,Ti], 

with c e G R 2 . By letting e —> in the previous relation we obtain a contradiction with the assumption 
that t — > 4>(t,£) is not constant for at least one £ G Xo, in fact to is any point of [0,T], and T = 
n T , n GN. □ 



5 Example 

In this section we consider the following illustrative example for system (2). 

21 = X 2 - Xl{x\ +X% - 1) , . 

x 2 = -x\ - x 2 {x\ + x\ - 1). 

It is easy to see that system (33) has the limit circle xq{6) = (sin 8, cos 0) with period To = 2ir. In 
order to verify the conditions of theorem 1 we use theorem 3. Thus, we consider 

*>« -(-£.')■ 

It is easy to verify that y(t) = e -2t (sint, cos t) satisfies the linearized system (3) corresponding to 
system (33). Therefore 

-1 _ „-2t / - COS t \ 



# = e 
y(t) T = e 



sint 



T _ -2t / COSt 



sint 
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Define the function / : [0, 2ir] — > (R) by the formula 

/(*)= / 



cos 9 



If £ — > 4>(t, x) is a 47r-periodic function, then condition (B2) of theorem 3 takes the form 

7 ( COST e -t TSinT Y\(T-e, X0 (T))dT,( Sin ^W0 

y \ Sill t e COS T / V cos 1/ 



s s-4tt+9 



for any s G [0,4vr], G [0,2vr]. 

As example of function <\> satisfying condition (34) we consider 



<t>(t,0 



£2 Ci 1 / £1 cos i — £2 sin i + a sin 



-£l 6 y V ^1 sin t + ^ 2 cos t 
with a > sufficiently small. In fact, for a = we have 



/" / cos r e sin r \ , . _ , . . , / sin ( 

/ • -2r <p[T -6,Xa(T))dT, 

J \ — sin t e cos r / 



cos ( 



/ cos r — sin r \ / cos r sin r 
\ e 2r sin r e 2r cos r / \ — sin r cos r 

(sin r cos(r — 6) — cost sin(r — 6) \ / sin0 
sin r sin(r — 0) + cost cos(t — 6) J 'I cos0 

= /(s,n«cos«) T f I "U/*"^ 



e 2r / ' I cos 



s-4tt+6» 

s+6» 



= 4vr sin 2 + cos 2 ^ e 2r (iT > 



s—4ir+9 



and so (34) holds for the function (35) for any a > sufficiently small. 
Let us now verify condition (B3) of theorem 3. For this, consider 



N(9) 

and observe that the homotopy 

N x (0) 



e 29 cos sin \ / sin 
-e~ 26l sin0 cos0 / I cos0 



e 2X9 cos sin \ / sin ( 
-e~ 2A6l sin0 cos0 / I cosl 



joining the vector fields Nq and N\ = N does not vanish whenever A G [0, 1]. Therefore 

sin(2- 
cos(2- 



deg(A, [0, 2tt]) = deg(A , [0, 2tt]) = deg ( { ] , [0, 2vr] 
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and so condition (B3) of theorem 3 is also satisfied. 

In conclusion, for < e < e 7 , where e 7 > is given as in (8), the following system 




x 2 - x\{x\ + x\ — 1) 
-Xl - x 2 (x? + x| — 1) 



2 ^ 1 oj H 

xi cos t — X2 sin i + a sin | \ 
xi sint + X2 cos t y ' 

where a > is sufficiently small, has a 47r-periodic solutions x £j i and x £j 2 which converge to the unitary 
circle from the outside and from the inside respectively as e — ► 0. 
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